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In this study, we propose new Higgs production mechanisms with multi-photon final states in the 
fermiophobic limit of the Two Higgs Doublet Model. The processes are: gg — > hh, gg — > Hh 
followed by H — * hh and gg — > Ah followed by A — * hZ. In the fermiophobic limit, gg — > hh 
and gg — > Ah — > hhZ would give rise to 47 signature while gg — > Hh — > hhh can give a 67 final 
state. We show that both the Fermilab Tevatron and CERN's Large Hadron Collider can probe 
a substantial slice of the parameter space in this fermiophobic scenario of the Two Higgs Doublet 
Model. If observed the above processes can give some information on the triple Higgs couplings 
involved. 



I. INTRODUCTION 



There are several extensions of the Standard Model (SM) with an enlarged scalar sector. Some of these 
extensions allow for Higgs with reduced or even no couplings to the fermions. They are referred to as 
fermiophobic Higgs scenarios in the literature [l[ . The DO collaboration has recently presented new results 
on fermiophobic Higgs searches @, 0| ■ In Q they have searched for a fermiophobic Higgs produced in 
association with a charged Higgs. The full process is pp — > hH ± — > hhW^ * — > 47 + X and was proposed 
in [1, 0, 0] . DO required at least three photons in the final state for maximizing signal efficiency. For each 
pair of values of (tan (3, m H ± ) a bound on the fermiophobic Higgs mass was set. In DO has performed 
a search for the inclusive production of di-photon final states via the Higgsstrahhmg and vector boson 
fusion processes: pp — > hV — > 77 + X and pp — ► VV — ► h — > 77 + X, respectively with a total integrated 
luminosity of 1.10 ± 0.07 /£> -1 . A lower nih bound of 100 GeV was obtained in a benchmark scenario 
that assumes hVV (V= W, Z) couplings to be exactly the same as in the SM and all fermion branching 
ratios to be exactly zero. 

All LEP collaborations have searched for a fermiophobic Higgs in the channel e + e~ — > h(— > 77) Z 
@j H, H 03 ■ The combination of all results [ll[ yielded the lower bound for the fermiophobic Higgs 
mass of 109.7 GeV, at 95% confidence level, which again is valid only in the above benchmark scenario. 
Searches in the e + e~ — > hA channel were also performed at LEP with lower bounds derived for nth + rriA 
(see references 0, Q for details). The new DO bound Q on is weaker but a larger region of the 
model's parameter space is covered. The channel qq' — > V* — ► h(— ► 77) V had already been used at the 
Tevatron by the CDF [HI and by the DO [U[ collaborations to set limits of 78.5 GeV and 82 GeV at 
95% confidence level, respectively, on the fermiophobic Higgs mass. It is expected [H, GJ, [lj] that all 
Tevatron bounds will improve once the data collected at 2/6 -1 luminosity is analyzed. 
There are however other ways of producing a four photon final state in a fermiophobic scenario. In this 
letter we will consider all fermiophobic Higgs production processes with at least three photons in the 
final state produced via gluon fusion. As shown in Q, a signal with at least three photons is very easy 
to extract at the Tevatron. The most relevant process for the analysis is gg — » hh — > 4 7 but we can also 
have gg — > hH — > hhh -^67 and gg — ► hA — > hhZ — > 4 7 + X. Wc will show that both the Tevatron and 
the Large Hadron Collider (LHC) can probe a substantial region of the parameter space. We will discuss 
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as well the complementarity between the different production modes. 

This paper is structured as follows: we will review the fermiophobic model in section II and then proceed 
to look at the all available theoretical and experimental bounds in section III. In section IV, we will discuss 
in detail the production process and in section V the signal. Analysis of the results and conclusions will 
be presented in section VI. 



II. THE FERMIOPHOBIC THDM 



To define our notation we start with a brief review of the two-Higgs doublet potential used here. The 
potential chosen is the most general, renormalizable, CP-conserving potential, invariant under SU(2) (g> 
U(l) that one can build with two complex Higgs doublets. It can be written as 

V(^, $ 2 ) = Axd^l 2 - vlf + A 2 (|$ 2 | 2 - v 2 ) 2 + A 3 ((|3>i| 2 - vl) + (|$ 2 | 2 - v 2 )) 2 + 

A 4 (|$i| 2 |$ 2 | 2 - |*i^ 2 | 2 ) + A 5 [5Re($+$ 2 ) - v lV2 ] 2 + A 6 [3to($+$ 2 )] 2 (1) 

where $1 and $ 2 have weak hypercharge Y=l, v\ and w 2 are respectively the vacuum expectation values of 
<3?i and < & 2 and the Xi are real-valued parameters. Note that this potential violates the discrete symmetry 
<f>; — > — $; only softly by the dimension two term A55Re($^<i> 2 ). The hard breaking terms (dimension 
four) of the discrete symmetry have been set to zero. As in all other THDM, we end up with two CP-even 
Higgs states usually denoted by h and H, one CP-odd state, A and two charged Higgs bosons, i? ± . The 
potential in eq. {1} has 8 parameters (including v\ and w 2 ). The combination v 2 = v\ + v\ is fixed as 
usual by the electroweak breaking scale through v 2 = (2v / 2G^ We are thus left with 7 independent 
parameters; namely (Xi)i=i,...,e, and tan/3 = v-ijv\. Equivalently, we can take instead 

TUh , ttih , m a , itlh± , tan/3 , a and A5. (2) 

as the 7 independent parameters. The angle (3 is the rotation angle from the group eigenstates to mass 
cigenstates in the CP-odd and charged sector. The angle a is the corresponding rotation angle for the 
CP-even sector. 

In a general THDM it is possible to couple just one doublet to all fermions by choosing an appropriate 
symmetry for both the fermions and the scalars. This model is known as THDM type I in the literature. 
Like in the SM, where just one doublet couples to all fermions, each scalar couples to the different 
fermions with the same coupling constant. However, unlike the SM, the couplings are now proportional 
to the rotation angles a and /3. For instance, the lightest CP-even Higgs couples to the fermions as 
cos a. I sin/3 gf¥f- By choosing cos a = 0, the lightest CP-even Higgs decouples from all fermions. It is 
usually referred to as a fermiophobic Higgs scalar This way the heavy CP-even scalar will acquire 
larger couplings to the fermions than the corresponding SM couplings. The remaining scalars are not 
affected by this choice as they do not couple proportionally to a. However, h can still decay to two 




FIG. 1: Feynman diagram of the largest contribution to h — » bb 



fermion pairs via h W*W{Z*Z) -»■ 2// or h -> W*W*{Z*Z*) 2 //. We will include these decays 
in our analysis. It is worth pointing out that these processes occur near the W(Z) threshold. Decays of 
h to two fermions can also be induced by scalar and gauge boson loops (see e.g. fig. [T]). In the THDM, 
the angle a has to be renormalized to render h —> f f finite. However, at a — ir/2, all one-loop decays 
h — » // are finite. Thus we can impose the following condition for the renormalization constant 8a: the 
renormalized one-loop decay width for h — > // is equal to the finite unrenormalizcd decay width. This 
condition is equivalent to setting [5a] a=7T / 2 — 0. In [131 we nave checked that this condition holds for all 
fermions. The only relevant one-loop decay is h — > bb due to a large contribution of the Feynman diagram 
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shown in fig. [T] to the total decay width j28[. Thus, on one hand, h is not completely fermiophobic at 
a = 7r/2, and on the other hand, all decays h — ► ff but h — > bb are almost zero even at one-loop level. 
Nevertheless it is possible, although hard, to find regions of the parameter space where h — > bb has a 
sizeable effect in the fermiophobic Higgs signature. With this in mind it is obvious that we are working 
with a version of the model with one less parameter than the general CP-conserving THDM. 



III. THEORETICAL AND EXPERIMENTAL BOUNDS 

In our analysis we take into account the following experimental and theoretical constraints. Note however 
that we will use the theoretical constraints as guides and explain in each case how they would affect the 
results given in the plots. We believe that if Nature did not provide a fermiophobic Higgs it should be 
disallowed ultimately by experiment. 

• Experimental bounds from LEP; the LEP collaborations have set bounds on the mass of a 
fermiophobic Higgs as a function of tan (3. The most stringent bound comes from the combination 
of all LEP experiments given in [ll|. The 95% CL limit on tan/3 is about 18 for run = 20 GeV. 
From nih = 20 GeV until rrih = 70 GeV, the bound oscillates about tan/3 = 10. For nth > 70 
GeV, tan/3 > 10 is already a conservative bound. In the plots tan/3 = 10 is used because much 
higher values of tan /3 would violate pcrturbativity constraints. Note however that as we will show, 
gg — > hh is independent of tan (J and for tan (J > 10 the dependence of all other relevant processes 
is negligible. 

There is another bound on the fermiophobic Higgs coming from hA production at LEP which 
constrains severely the value of tan (3 especially for small h masses. If however we take the A mass 
to be above 150 GeV, the hA production is no longer a constraint for all values of tan/3. 
Finally, if tan/3 is large, the ZZh coupling is suppressed while the non-fcrmiophobic CP-even 
Higgs H will couple to the Z bosons with almost the SM strength. Therefore, in the plots shown, 
the minimum value of the H mass is 100 GeV and most cases presented are for an H mass above 
120 GeV. 

• As already stated in the introduction, the DO and CDF experiments recently reported searches 
for a fermiophobic Higgs in two different channels. DO searched 0] for the inclusive production of 
di-photon final states via the Higgsstrahlung and vector boson fusion. The bound with the model 
benchmark described in the introduction is weaker than the LEP bound but it spans a larger region 
of the parameter space. The search Q in the pp — > hH ± — > hhW^ * — > 4j + X channel sets a 
bound on the fermiophobic Higgs mass for each pair of values of (tan /3, m H ± ) . We will take all 
these bounds into consideration in our analysis. 

• The extra contributions to the dp parameter from the Higgs scalars [HI should not exceed the 
current limits from precision measurements [l9j : \Sp\ < 10~ 3 . Such an extra contribution to Sp 
vanishes in the limit m H ± = ttia- To ensure that Sp will be within the allowed range whenever 
possible we allow only a small splitting between m H ± and ttla- 

• Recently, it has been shown in Ref. [2(| that for THDM models of the type II, data on B — > X s ^ 
imposes a lower limit of mu± > 290 GeV. In THDM type I, there is no such constraint on the 
charged Higgs mass. Therefore, in our numerical analysis which is valid for THDM type I we will 
ignore the limit on the charged Higgs. 

• The scalar sector can also be constrained using perturbativity constraints on A^ [2l|, [2^] ■ In the 
present study we will not impose those constraints in order to quantify the optimal cross sections 
and scan over all parameters space. In fact, in the fermiophobic limit, the process 

99 -> hh 

for example depends only on m/j, rati and A5. As we will explain later, the tan/3 dependence in 
gg — > hh drops out. As a result, for a given mh, tuh and A5, one can tune tan/3, tua and mj?± in 
order to satisfy perturbativity constraints. 

• From the requirement of perturbativity for the top and bottom Yukawa couplings [23|, tan/3 is 
constrained to lie in the range 0.3 < tan/3 < 100. But it turns out that from pcrturbativity 
argument on A,, moderate values of tan/3 less than about 10 are preferred. 
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In order to respect perturbativity constraints we will use moderate values for M^, Mh, Ma and M H ±. 



IV. PRODUCTION CROSS SECTIONS 



The process pp(p) — ► hh has both tree level contributions mediated by Higgs exchange from qq — > 
H*,h* — ► /i/i and one loop contributions from gluon fusion gg — ► /i/i. The tree level contribution is 
proportional to the quark masses which will be neglected for the Tcvatron energies. We have also checked 
that even in the large tan (3 limit the production process qq — > hh is also negligible for the LHC when 
compared to gg — > hh. The process gg — > hh occurs only at the one-loop level. As we will see, even 
if loop suppressed, this process can still be enhanced by the strong QCD coupling as well as by the 
heavy Higgs H resonant effect when it can decay to two light CP even h scalars. There are two types 
of diagrams that participate in the process gg — > hh. The box diagram, as shown in fig [2] represents 
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FIG. 2: Box contribution to gg — > hX with X = H,h, A and / is a fermion. Just one topology is shown. 



just a generic contribution to the process. We have also included all other quarks and because both 
the initial and the final state have identical particles, there is a total of six diagrams for each flavor. 
The second type are the vertex diagrams which again are six for each flavor. In fig. [3] we show just the 
representative diagrams with a generic fermion in the loop. In the fermiophobic limit, the top loop is 
always the dominant contribution. The other two processes, pp(p) — > H h and pp(p) — > Ah have a similar 




FIG. 3: Triangle contribution to gg — » hX with X = H,h, A. 



structure (vertex and box) to pp — > hh, except for one additional contribution to the vertex diagrams in 
PP(P) ~~ * Ah which is the Z boson and the Goldstone boson s-channcl exchange. 

However, in the fermiophobic limit, all box contributions vanish and the same is true for all vertex 
contribution with s-channcl Higgs (h) exchange due to the fact that the fermiophobic Higgs coupling 
to fermions is zero. Since the box contributions drop out in the fermiophobic limit, the two processes 
PP{P) ~* hh, Hh are directly proportional to the pure scalar couplings Hhh and HHh. The third process 
PP(P) —* Ah is sensitive both to pure scalar couplings AAh and AGh as well as to the gauge coupling 
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ZAh. Hereafter we list the pure scalar coupling in the fcrmiophobic limit (FL) needed for our study: 



K Hhh 



\FL 
^HHh 



A AAh 



\ FL 
A AGh 



\FL 



eX$v 2 sin/3 tan/3 
oc A 5 - 



4m w s w + tan 2 



eX^v 2 cos/3 



oc A 5 



4m w s w + tan 2 /3 



Amwsw cos /3 



2 sin 2 /3 m 2 + X 5 v 2 - 4 cos 2 /3 



e sin/3 
e 



7714 — to. 



2mvi/svi/ sin 2/3 



2 sin 3 /3to 2 - X 5 v 2 sin /3 - 2 sin 2/3 cos ^m 2 ^ 



(3) 
(4) 
(5) 
(6) 
(7) 



Note that the Hhh and HHh couplings are directly proportional to A5 while the AAh coupling depend 
both on As, mj as well as on m/,. It is clear that in the case of exact discrete Zi symmetry (A5 = 0), 
both Hhh and HHh would vanish. In the fermiophobic limit, all fermions couple to each scalar with the 
same strength. The ffH coupling is proportional to 1/ sin/3. Therefore from eq. [3] we conclude that the 
(3 angle dependence cancels out in the cross section, that is, a g g—>hh depends only on mh, rnn and A5. 
The ffA coupling is proportional to cos /?/ sin (3 which means that in the large tan/3 limit the /3 angle 
dependence is also very mild except for very large h and/or A masses. Those cases will not be included 
in our study. Finally, note that not only <J gg —>hh vanishes in the limit A5 = but also <j gg ->Ah, in the high 
tan (3 limit becomes negligible when A5 = due to to the smallness of the values of the masses involved. 
The one-loop amplitudes were generated and calculated with the packages FeynArts [24[ and FormCalc 
[25l ]. The scalar integrals were evaluated with LoopTools [26j |. 



1. Numerical results 



In this section, we present our numerical results. As stated earlier, to avoid the LEP bounds on the 
fermiophobic Higgs we will fix tan (3 to be of the order 10 for the processes gg — » Ah and gg — > Hh. This 
way we suppress the ZZh coupling while keeping perturbativity bounds on A; within the allowed range. 
The first consequence of this choice of tan/3 « 10 is that the coupling Hhh is enhanced (sin j3 « 1) while 
HHh is suppressed (see eqs. [3l [4]). As we have discussed, in the fermiophobic limit, the process pp — ► hh 
has only vertex contribution through s-channel heavy Higgs (H) exchange. Therefore, the cross section 
for pp{p) — > hh depend both on 777 h as well as on Hhh coupling which is proportional to A5. There are 
two sources of enhancements for pp — > hh. The first one is to take Hhh (or equivalcntly A5) as large as 
possible. The second one is when hh production is resonant, that is, mu w 2to^. 

We first discuss pp(p) — > hh production. In Fig. 2] (left) we illustrate the cross section of pp — > hh as a 
function of the fermiophobic Higgs mass mh for two representative values of Ag = 4ir and 871- . The other 
heavy CP even mass is taken to be to# = 2mh such that the resonant channel H — > hh is open. The 
other parameters are: 777,4 = tuh± = 300 GcV. On the H resonance, the cross section is enhanced and 
can reach a few hundreds of picobarn for a very light fcrmiophobic Higgs 777^ = 50 GeV. Even for mh of 
the order of 100 GcV, the cross section is still larger than 0.1 pb which would give thousands of produced 
events for the planned Tevatron luminosity of 10/6 -1 . As can be seen from the analytic expression of 
the coupling Hhh, the large A5 is the large is the coupling Hhh and so is the pp — > hh cross section. 
The observed kinks at mj? = 160 and 182 GeV in the left plot are due to the opening of the decay 
channels H — > WW and H — > ZZ. On the right-hand side of Fig. [4] we show the same cross section, 
the only difference being the H mass which is now fixed at 120 GeV. It is clear that even away from 
the H resonance a significant set of fermiophobic Higgs masses and A5 values can still be probed at the 
Tevatron. The behavior with A5 changes at threshold due to the H width effect on the cross section. 
Above threshold the H — > hh channel is closed, the H width is then very small. Therefore in that region 
the cross section is just proportional to A 2 . Below threshold, the H — > hh channel is open and the width 
in the H propagator starts to play a role which makes the dependence with A5 no longer trivial. 
In Fig. [S]we show the same plots as in Fig. 0] but for the LHC. As expected the plots are re-scaled by 
more than one order of magnitude. For the values shown, all masses between 10 and 100 GeV can be 
probed for most values of A5. Even for A5 = ir/2 the smallest cross section value is of the order of 10 fb. 
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FIG. 4: a(gg — > hh) for the Tevatron in units of pb as a function of the fermiophobic Higgs mass with As = 4ir 
and 8n. On the left the mass of the heavier CP-even Higgs boson is twice the fermiophobic Higgs mass, i.e., on 
threshold, and on the right it is fixed (120 GeV). 
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FIG. 5: u(gg — * hh) for the LHC in units of pb as a function of the fermiophobic Higgs mass with As = tv/2 , 4ir 
and 8n. On the left the mass of the heavier CP-even Higgs boson is twice the fermiophobic Higgs mass, i.e., on 
the threshold, and on the right it is fixed (120 GeV). 



In Fig. [6] we show the cross section of pp — > hh as a function of the resonant Higgs mass mn for A5 = 87r 
and for a fermiophobic Higgs mt = 60 GeV, for the Tevatron (left) and for the LHC (right). As one can 
see, away from the resonance mu as 2nih, the cross section is of the order a few fb for the Tevatron and 
a few pb for the LHC. Once we cross the resonance, one can see a spectacular pick for run = 120 GeV 
which is due to the opening of the decay channel H — > hh. The peak is very sharp because the total 
width of H is very narrow for ma < 2m h, less than 10~ 2 GeV. Once the decay channel H — > hh is open 
for mn <; 2mh the decay width of H increase suddenly to more than 100 GeV. This is manifestly seen 
in the plot by a dramatic decrease of the cross section from few hundred pb to 10~ 2 pb. 

Finally we present the pp(p) —* Hh and pp{p) — > Ah reactions. It is clear from both plots presented in 
Fig. [7] that the cross section for Hh production is negligible for most of the parameter space. On the 
contrary, the cross section for hA production can be very large and still within the Tevatron reach. The 
cross section gg — > Ah can be several orders of magnitude larger than the corresponding qq — > Z* — > Ah 
(see ref. We will show in the next section that taking into account the pseudo-scalar branching 

ratios, the decay A — > Zh can be the dominant one, and so the channel gg —> Ah —> Zhh — > Z47 is still 
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FIG. 6: a(gg — > hh) for the Tevatron (left) and for the LHC (right) in units of pb as a function of the heavier 
CP-even Higgs mass with A5 = 8n and a fermiophobic Higgs mass of 60 GeV. 
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FIG. 7: <y(gg — * hA) and o(gg — » /i-ff) for the Tevatron (left) and for the LHC (right) in units of pb as a function 
of the fermiophobic Higgs mass with tan (3 = 10, As = 47r, 87r, Mh — 150 GeV and Ma = 150 GeV. 



worth exploring at the Tevatron. Finally, we have checked that the contribution from qq — > hh, Hh Ah 
was negligible. 



V. HIGGS SIGNATURE 



Having established that gg — * hh and gg — > ^4ft. are worth studying both at the Tevatron and at the LHC 
we now turn to the experimental signatures for the fermiophobic Higgs and for the pseudo-scalar in the 
parameter space under study. In the fermiophobic limit there is a dramatic change in the fermiophobic 
Higgs signatures. For smaller fermiophobic Higgs masses, the main decay is to two photons (through W 
and charged Higgs loops) until the WW channel starts to dominate. The crossing point of the branching 
ratios depends on the parameters of the scalar potential which enters the game through the charged 
Higgs contribution to h — * 77. Those parameters are mainly the charged Higgs mass, A5, tan/3 and 
the fermiophobic Higgs mass (see eq.©). Recently a detailed study of h — ► 77 appeared in [27|. Due 
to all experimental and theoretical constraints, for a fermiophobic Higgs with mass between 10 and 100 
GeV, Br(h — > 77) w 100% except in a tiny neighborhood of A5 = 0. In this neighborhood, all cross 
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sections are extremely small. Therefore there will always be a tiny region around A5 = that will not be 
probed with the processes proposed here. For illustration, we show in fig. [8] (left) the fcrmiophobic Higgs 
branching ratio as a function of A5 . We have checked that the larger tan /3 is the higher the value of the 
77 branching ratio. Br(h — > 77) decreases with the fcrmiophobic Higgs mass, but for rrih = 100 GeV 
the plot looks almost the same. Regarding the pseudo-scalar decays, as one can see in fig. [5] (right), one 
can have Br (A — > hZ) ss 100% if the decay A — > H^W^ is kinematically forbidden. We have checked 
that changes in tan/3, A5 and the A mass produce negligible changes in the branching ratio provided the 
charged Higgs channel A — > H^W^ is kept closed. Therefore, the 4 7 final state is by far the dominant 
one. 



VI. DISCUSSION AND CONCLUSIONS 



In this work we have shown that there are alternative channels to search for fermiophobic Higgs with a 
multi-photon signature. A vast region of the parameter space of the fermiophobic THDM can be probed 
at the Tevatron and the analysis can easily be extended for the LHC. In the fcrmiophobic limit, the 



angle (3 is already very constrained. LEP [llj has set a limit of tan/3 > 10 for almost all values of the 
fermiophobic Higgs masses up to 100 GeV. For some masses the bound is even stronger. On the other 
hand theoretical constraints tell us that these values cannot be too high. This implies that gg — -> hh and 
gg — > hA will be large while gg — > hH will be negligible. Note however that the gg — ► hh cross section 
does not depend on tan/3 and that the gg — > hA dependence on tan/3 for values above 10 is negligible. 
Regarding the A5 dependence, we have shown that there is a "tiny to small" region around A5 = that 
can not be probed. This is especially true for hh productions but gg — > Ah decreases with A5 as well. 
The process pp — > hH ± -> hhW ± * -> 47 + X proposed in 0, i] and studied in Q is complementary to 
the processes we propose in this work. We probe the region of small H and/or A masses while Q probes 
the small charged Higgs mass region. The advantage in our case is that our study is independent of tan /? 
while their process does not depend on A5 . Finally if all Higgs scalars besides the fcrmiophobic Higgs are 
very heavy, only the two photon search can exclude a fermiophobic Higgs. 

As expected and as it was shown in [2| the background for a 3 or 4 photon final state is easy to control. 
In the case of the 3 photon final state the main background contribution comes from the direct tri-photon 
production (see for details). 
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